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I .   INTRODUCTION 

The  earth's  atmosphere  has  long  been  of  interest  to 
science.  From  explaining  phenomenon  such  as  the  aurora 
borealis  or  northern  lights  to  the  study  of  weather  and 
climate,  knowledge  of  the  atmosphere  has  many  interesting  and 
useful  applications.  In  addition  to  the  scientific  community, 
the  military  considers  an  understanding  of  the  atmosphere  to 
be  of  extreme  importance  because  of  its  profound  effect  on 
several  critical  systems.  Specifically,  the  characteristics 
of  the  region  of  the  upper  atmosphere  known  as  the  ionosphere 
are  of  particular  significance.  Among  the  features  of 
interest,  the  electron  density  is  of  principal  concern  due  to 
its  role  in  high  frequency  (HF)  radio  communications, 
satellite  communications  and  control,  and  over-the-horizon 
radar  (OTHR).  Another  important  aspect  is  the  naturally 
occurring  luminescence  of  the  atmosphere  called  airglow.  This 
is  significant  because  it  offers  a  way  to  determine  electron 
densities.  Furthermore,  airglow  is  the  source  of  a  background 
noise  signal  that  affects  infrared  detection  and  tracking  of 
ships,  aircraft  and  ballistic  missile  launches  by  space 
systems.  The  electron  density  and  airglow  are  not  constant 
but  vary  due  to  several  factors  such  as  geographic  location, 
time  of  day  and  solar  activity. 


Currently,  electron  densities  in  the  ionosphere  are 
determined  by  actively  sounding  the  atmosphere  with  radio 
waves.  This  is  accomplished  in  the  following  way.  A  radio 
wave  of  a  certain  frequency  is  transmitted.  This  signal  is 
reflected  by  the  ionosphere  when  it  encounters  a  unique 
electron  density  which  corresponds  to  the  transmitted 
frequency.  This  frequency  in  called  the  critical  frequency. 
The  reflected  signal  is  detected  by  a  receiver  and  the  time  it 
takes  to  make  the  round  trip  is  recorded.  By  multiplying  half 
this  time  by  the  velocity  of  waves  in  free  space,  one  can 
determine  what  is  known  as  the  "virtual"  height.  This  is  not 
the  actual  height  at  which  the  reflection  occurred.  The 
calculated  virtual  height  is  actually  greater  than  the  true 
height  of  reflection.  This  is  because  the  wave  group  velocity 
in  the  ionosphere  is  less  than  the  free-space  group  velocity. 
Techniques  exist  for  determining  the  true  reflection  height 
from  the  virtual  height.  By  varying  the  frequency,  an 
al  titude-versus-density  profile  can  be  obtained.  These 
soundings  are  conducted  by  instruments  called  ionosondes  and 
while  they  are  accurate,  they  have  significant  disadvantages. 
First,  these  are  fixed  sites  with  limited  coverage  and  may  not 
be  located  in  a  particular  area  of  interest.  Second,  radio 
sounding  is  only  capable  of  providing  the  lower  half  of  the 
electron  density  profile  from  the  ground.  This  is  because 
although  the  same  electron  density  may  occur  at  different 
altitudes,  the  radio  wave  will  always  reflect  off  the  first 


appropriate  density  it  encounters  i.e.,  the  one  which  is  at 
the  lower  altitude.  Finally,  ionosondes  do  not  provide  a 
measurement  of  the  airglow  intensity  or  densities  of  other 
constituents  of  the  ionosphere. 

As  an  alternative,  one  might  suggest  putting  an  ionosonde 
on  a  satellite.  This  might  correct  the  coverage  problem  but 
would  then  only  provide  the  upper  half  of  the  density  profile 
and  would  still  not  give  information  about  the  airglow  or  the 
additional  densities  of  interest.  Furthermore,  due  to  its 
active  nature  (i.e.,  it  transmits  a  signal),  an  ionosonde 
would  necessarily  be  large  with  a  large  power  requirement  and 
thus  not  well  suited  for  satellite  deployment. 

An  alternative  solution  is  to  place  a  spectrograph  on  a 
satellite  which  could  record  emissions  from  the  airglow. 
These  emissions  could  then  be  used  to  infer  the  densities  of 
neutral  and  charged  species  from  which  the  electron  density 
can  be  determined.  Because  it  is  a  passive  instrument,  a 
spectrograph  is  less  complex,  requires  less  power  and  is 
smaller  than  an  active  system  making  it  ideal  for  satellite 
deployment.  One  such  instrument  is  the  Naval  Postgraduate 
School's  middle  ultraviolet  spectrograph  (MUSTANG).  The 
MUSTANG  was  part  of  a  scientific  payload  in  two  successful 
rocket-borne  experiments:  one  conducted  in  March  1990  and  the 
other  in  February  1992.  This  thesis  analyzes  the  data 
collected  from  the  1990  flight. 


A.  THESIS  OBJECTIVES 

It  is  the  purpose  of  this  thesis  to  present  a  method  of 
determining  the  densities  of  0,  02  and  Nj  by  comparing  data 
with  the  model  calculations.  Additionally,  the  procedure  for 
determining  the  densities  of  the  charged  species  and  electrons 
will  be  addressed. 

B.  THESIS  OUTLINE 

This  thesis  is  divided  into  five  chapters.  Following  the 
introduction,  the  second  chapter  gives  background  information 
necessary  to  understand  the  experiment  and  the  analysis. 
Included  in  this  chapter  is  a  discussion  of  the  earth's 
atmosphere,  the  airglow,  and  the  atmospheric  models.  The 
third  chapter  is  a  brief  description  of  the  experiment 
including  an  explanation  of  the  MUSTANG  instrument's 
construction  and  how  it  works.  Chapter  IV  contains  the 
analysis  of  the  data  and  the  results.  The  final  chapter 
discusses  conclusions  and  suggestions  for  further  research. 


I I .    BACKGROUND 

Before  one  can  understand  the  analysis  in  Chapter  IV,  it 
is  important  to  be  familiar  with  the  concepts  and  models  on 
which  this  analysis  is  based.  This  chapter  provides  the 
necessary  background  knowledge  required  to  follow  the 
analysis . 

A.   THE  EARTH'S  ATMOSPHERE 

The  earth's  atmosphere  extends  from  the  surface  out  to 
approximately  1000  km.  Figure  2-1  shows  a  temperature  versus 
altitude  profile  of  the  atmosphere.  One  can  see  that  the 
atmosphere  is  made  up  of  several  temperature  gradients.  These 
gradients  can  be  thought  of  as  layers  describing  various 
temperature  regions.  The  layer  at  the  earth's^  surface  is 
called  the  troposphere.  Its  temperature  is  influenced 
greatest  by  the  surface  temperature  of  the  earth.  This  region 
is  characterized  by  a  negative  temperature  gradient.  The 
temperature  continues  to  decrease  until  an  altitude  of 
approximately  10  km.  At  this  point,  called  the  tropopause, 
the  temperature  starts  to  increase  with  altitude. 

This  next  region  is  called  the  stratosphere.  .The  rise  in 
temperature  in  the  stratosphere  is  due  to  heating  caused  by 
the  earth's  ozone  layer  absorbing  ultraviolet  light.  This 
region  has  received  much  attention  in  recent  years'because  of 
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Figure  2-1  Temperature  versus  altitude  profile  of  the  earth's 

atmosphere . 


the  possible  ozone-depleting  effects  of  chl orof 1 orocarbons 
(CFC's)  in  the  atmosphere.  These  man-made  chemicals,  used  as 
solvents  and  refrigerants,  and  in  many  industrial  processes, 
are  known  to  act  as  a  catalyst  in  a  reaction  that  destroys 
ozone.  Without  ozone,  the  earth  would  not  be  protected  from 
the  deadly  effects  of  ultraviolet  radiation  and  most  life 
forms  could  not  survive.  The  ozone  heating  peaks  at  the 
stratopause  at  an  altitude  of  around  50  km.  The  next  region 
is  another  negative  temperature  gradient  layer  called  the 
mesosphere.  Temperature  decreases  until  the  mesopause  is 
reached  at  an  altitude  of  about  86  km.  Above  the  mesopause  is 
the  thermosphere  where  temperature  increases  until  it  becomes 
isothermal  with  increasing  altitude.  This  isothermal 
temperature  is  referred  to  as  T,  and  varies  between  1000  and 
1800K.  This  isothermal  region  is  called  the  exosphere.  The 
exosphere  is  characterized  by  extremely  low  densities  and  is 
the  beginning  of  what  one  might  call  the  "vacuum"  of  space. 
The  region  of  interest  for  this  thesis  is  the  thermosphere. 

The  major  constituents  of  the  thermosphere  are  molecular 
nitrogen,  molecular  oxygen  and  atomic  oxygen.  Their  density 
profiles  and  the  temperature  profile  in  the  thermosphere  make- 
up what  is  referred  to  as  the  neutral  atmosphere.  Radiation 
from  the  sun  ionizes  these  species  in  a  process  called 
photoionization .  The  total  ionization  rate  for  a  species  at 
an  altitude  z,  Q(z),  is  determined  by  the  following  equation: 


Qa(z)=na(z)   f  o5(X)  Jx(oo)e-T(A^)CT(6)dX    (2-1) 


LCj= 


where : 

n* (z)  is  the  volume  density  of  species  a, 

of  (X)  is  the  ionization  cross  section  for  species  a, 
Jx(»)  is  the  solar  flux  at  the  top  of  the  atmosphere, 

t(z)  is    the  optical  depth  at  altitude  z, 

CH(Q)      is      the     Chapman    function. 

X^     is    the      threshold   wavelength 

The  optical  depth  term  reduces  the  total  ionization  rate  at  an 
altitude  by  accounting  for  the  absorption  in  the  atmospheric 
column  above.  The  Chapman  function  accounts  for  the 
increasing  column  density  as  the  sun  moves  from  the  zenith  to 
the  horizon. 

In  addition  to  creating  ions  and  free  electrons,  photons 
from  the  sun  can  create  excited  species  and  energetic 
electrons.  These  energetic  electrons  are  called 
photoelectrons .  Photoelectrons  can  cause  additional 
ionizations  and  excitations.   Also,  molecules  can  dissociate 


as  a  result  of  photoabsorption  or  photoel ectron  impact. 
Chemical  reactions  also  take  place  in  the  thermosphere  which 
create  other  minor  neutral  and  ion  species.  These  species  can 
also  be  ionized  and  excited.  All  of  the  above  processes  are 
called  productions.  Productions  of  one  species  result  in  the 
loss  of  others.  Losses  can  also  occur  when  dissociated 
molecules  recombine  or  when  excited  particles  relax  back  to 
lower  energy  states.  The  region  formed  by  these  reactions  is 
called  the  ionosphere. 

The  ionosphere  is  comprised  of  four  regions  D,  E,  Fl  and 
F2  .  These  regions  are  based  on  electron  density.  The  density 
of  each  region  changes  with  the  time  of  day,  season  and  solar 
activity.  The  D-region  extends  from  an  altitude  of  50  to  90 
km,  the  E-region  from  90-140  km  and  the  F-region  from  140- 
1000  km.  The  species  in  the  ionosphere  obey  the  one 
dimensional  mass  continuity  equation  (MCE)  which  is  the 
subject  of  the  following  section. 

B.   THE  MASS  CONTINUITY  EQUATION 

The  densities  of  all  of  the  constituents  in  the 
ionosphere  at  a  given  altitude  obey  the  following  mass 
continuity  equation  (MCE) : 

§£W=P-nL-**^L.  (2-2) 

at  ot 


The  term  on  the  left-hand  side  of  the  equation  is  the  change 
of  density  of  a  species  at  a  given  altitude  with  time.  The 
first  term  on  the  right-hand  side  of  the  equation  is  the 
production  processes  discussed  above.  The  second  term  on  the 
right  is  the  losses  discussed  above.  The  final  term  on  the 
right  is  the  flux  term  which  accounts  for  diffusion  of 
particles  from  outside  the  region.  For  most  species,  the 
entire  MCE  is  required  to  calculate  the  density  of  a  species 
at  a  given  altitude.  Some  species  are  created  and  destroyed 
so  rapidly  that  there  is  no  time  for  them  to  diffuse  to 
another  altitude  interval.  As  a  result,  the  flux  term  can  be 
ignored.  In  addition,  for  some  species  production  and  loss 
are  in  quasi-equil  ibrium,  so  the  left-hand  side  of  the 
equation  is  zero.  Under  these  conditions,  one  can  solve  the 
MCE  for  the  volume  density: 

n(z)*£.  (2-3) 

Li 

This  situation  is  called  photochemical  equilibrium  (PCE).  In 
the  discussion  below  the  atmospheric  emission  model  uses  PCE 
in  the  procedure  that  calculates  volume  emission  rates.  It 
calculates  the  volume  emission  rate  from  the  volume  excitation 
rate.  This  is  done  in  the  following  way.  If  we  assume  that 
the  entire  density  of  an  excited  species  relaxes  giving  off  a 
photon,  then  the  volume  emission  rate  is  the  density  times  the 
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Einstein  coefficient,  where  the  Einstein  coefficient  A  is  one 
over  the  lifetime  of  the  species.   The  equation  is: 

VE=nA  (2-3a) 

Now  substituting  equation  (2-3)  for  the  density  and  setting  L 
equal  to  A  (because  the  losses  are  due  to  relaxation)  we  are 
left  with  the  volume  emission  rate  equal  to  the  production. 
Since  each  species  produced  is  assumed  to  relax  and  give  off 
a  photon,  then  the  volume  excitation  rate  is  equal  to  the 
volume  emission  rate.  There  are  some  other  losses  such  as 
quenching  which  must  also  be  included  but  they  do  not  change 
the  above  argument  because  they  are  simply  dimensionl ess  terms 
that  reduce  the  magnitude  of  the  production.  Now  that  the 
atmosphere  and  the  MCE  have  been  discussed,  another  important 
topic,  airglow,  must  be  addressed. 

C.   THE  AIRGLOW 

As  stated  in  the  introduction,  airglow  is  the  naturally 
occurring  luminescence  in  the  atmosphere.  Airglow  is  always 
present,  day  and  night.  In  fact,  the  airglow  for  each  of 
these  time  periods  has  its  own  name.  Not  surprisingly,  they 
are  the  dayglow  and  the  nightglow. 

One  source  of  the  airglow  is  the  ionosphere.  When  the 
excited  atoms  and  molecules  of  the  ionosphere  relax  to  lower 
energy  states,  they  give  off  energy  in  the  form  of "-  a  photon 
according  to  the  equation: 
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tf=hv  =  ^£.  (2-4) 


The  different  species  have  several  different  energy  levels 
within  a  state.  These  slightly  different  energy  differences 
result  in  a  range  of  wavelengths.  This  range  of  wavelengths 
is  called  an  emission  band.  When  an  intensity-versus- 
wavelength  profile  is  made,  the  result  is  a  spectrum.  In  a 
simple  case,  only  a  single  atomic  or  molecular  emission  is 
present  for  a  given  wavelength  range.  If  this  were  the  case, 
then  the  total  intensity  of  the  emission  band  can  be 
determined  by  integrating  over  the  wavelength  range  of  the 
spectrum.  Figure  2-2  shows  a  schematic  example  where  one 
emission  band  produces  the  entire  spectrum.  Unfortunately, 
nature  is  much  more  complicated.  A  spectrum  is  the  sum  of  the 
spectra  of  several  processes.  These  processes  produce 
emission  bands  that  overlap  or  contaminate  one  another. 
Figure  2-3  shows  the  realistic  case  of  a  spectrum  resulting 
from  the  sum  of  several  contaminating  emission  bands.  A 
technique  called  synthetic  spectral  fitting  can  be  use  to 
determine  the  intensity  for  the  individual  emission  bands. 
This  is  described  in  Chapter  III. 
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Figure  2-2   Simple  schematic  spectrum 
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The  1990  MUSTANG  data  was  collected  during  the  day, 
therefore  the  emission  bands  used  to  analyze  the  data  are 
dayglow  emission  bands.  One  of  these  emissions  was  the  N2 
(0,5)  Vegard-Kaplan  (VK)  band.  The  transition  for  the  VK  band 
system  is  AT,,  -  X"Z  ..  This  is  a  forbidden  transition  because 
excitation  from  the  ground  state,  X,  to  the  A  state  is 
forbidden  by  the  selection  rules  of  quantum  mechanics.  This 
does  not  mean  that  this  transition  cannot  occur.  It  means 
that  the  A  state  cannot  be  populated  by  photoexcitation .  The 
A  state  is  populated  in  two  ways:  directly  by  photoel ectron 
impact  and  indirectly  by  cascading  from  the  C,  B  and  W  states 
(Meier,  1991).  Another  emission  band  observed  in  the  data  is 
N2  Lyman-Birge-Hopf  ield  (alllg  -  ^2  )  .  This  is  another 
forbidden  transition.  The  a  state  is  populated  by 
photoel ectron  impact  (Meier,  1991).  Unlike  the  VK  bands,  this 
is  the  only  production  source  for  LBH.  The  final  emission 
observed  in  the  data  was  the  2972A  line  of  atomic  oxygen. 
This  line  is  designated  as  OI2972A.  This  is  a  :S  -  3P 
transition  and  it  too  is  forbidden.  0(  S)  has  different 
production  and  loss  reactions  for  low-altitude  (80-115  km)  and 
high-altitude  (above  115  km).  In  the  low-altitude  regime,  the 
source  of  0(*S)  is  the  following  reaction: 

0  +  0  +  M-02,+M  (2-5a) 

02*  +  M  -  02  +  0(!S)  (2-5b) 

where  M  is  taken  to  be  either  N2 ,  02 ,  or  0  so  that  the  density 
[M]  is  equal  to  the  sum  of  the  densities  [N2] ,  [02],  and  [0]  . 
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The  sink  for  this  altitude  regime  in  addition  to  radiative 
relaxation  is: 

O^S)  +02(a1Ag)  ^0+02.  (2-6) 

The  source  reaction  has  a  rate  coefficient  of  4.7xlO"33  (300/T)2 
cm  s  .  Because  the  proportion  of  Oj  in  the  a  state  was 
unknown,  the  chemical  loss  term  (Equation  2-6),  is  not 
included  in  the  model  calculations.  This  does  not  seriously 
affect  the  analysis  because  the  low-altitude  regime  is  below 
the  region  of  interest.  The  source  reaction  also  has  little 
effect  in  the  region  of  interest  but  was  included  for 
completeness.  In  the  high-altitude  regime  the  source 
reactions  are: 

0+e*-0(15)+e  (2-7a) 

02*+e-0+0(15)  (2-7b) 

o;+N^NO^O(xS)  (2-7c) 

Equation  (2-7b)  has  a  rate  coefficient  of  2.1  x  10   cm  s  l   for 

a  0(1S)  yield  of  10%  (13=0.1).   Equation  (2-7c)  has  a  rate 

-11    ^  -l 
coefficient  of  5  x  10   cm  s  .   The  sinks  are  the  following: 

0(1S)^0(1D)+h\5511  (2-8a) 

0(15)-0(3P)+i2V2972.  (2-8b) 

The  Einstein  coefficients  for  these  two  sinks  are  Ac 507= 1.0 6  s"* 
and   A2g72  =  0-045   s"  .      The   above   reactions   and   their 
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corresponding  rate  coefficients  were  taken  from  the  review  by 
Torr  and  Torr  (1982).  The  next  section  describes  the  various 
models  used  in  this  thesis  including  how  to  use  them,  what 
they  produce,  and  how  they  use  the  above  information. 

D.   THE  MODELS 

This  thesis  uses  three  main  computer  models  in  its 
analysis.  These  models  are  the  MSIS-83  model  atmosphere 
(Hedin,  1983),  MODATM ,  and  the  atmospheric  emission  model 
which  was  developed  for  this  thesis.  Each  will  be  discussed 
below  beginning  with  the  MSIS-83  model. 

The  MSIS-83  model,  developed  by  Hedin  (1983),  computes  a 
neutral  atmosphere  for  a  given  set  of  conditions.  When 
running  this  model ,  the  user  is  prompted  to  enter  several 
inputs.  The  first  is  IYD  which  is  the  year  and  .date  in  the 
following  format:  year  number  and  Julian  date  i.e.  YYDDD. 
Next  is  GLAT  which  is  the  latitude  in  degrees.  Then  the  user 
is  prompted  to  enter  GLONG  which  is  the  longitude  in  degrees 
from  -180  to  180.  After  the  longitude,  one  enters  a  number 
for  F107A.  This  is  the  three  month  average  of  the  F10.7  cm 
flux.  The  F10.7  cm  flux  is  a  radio  signal  produced  by  the  sun 
which  can  be  detected  on  the  ground.  It  gives  an  indication 
of  the  sun's  activity.  Next,  the  user  is  prompted  to  enter 
F107  which  is  the  F10.7  flux  for  the  day.  The  next,  input  is 
AP  which  is  the  Ap  index.   This  is  a  measure  of  magnetic 
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activity.  After  AP ,  one  is  prompted  to  enter  XLT  which  is  the 
local  time  in  hours  on  a  24  hour  clock.  The  last  input  is 
FRACXR  which  is  the  x-ray  scaling  factor  and  is  another 
measure  of  solar  activity.  The  values  for  F107,  F107A,  AP  and 
FRACXR  are  available  in  NOAA  publications  of  solar-terrestrial 
data.  Once  these  parameters  have  been  entered,  MSIS-83 
computes  the  neutral  atmosphere.  The  neutral  atmosphere  is 
produced  in  2  km  increments  from  86  km  to  450  km  .  The 
densities  of  0,  C>2  ,  and  Nj  as  well  as  the  temperature  profile 
are  stored  in  a  file  named  YYDDDMSIS.DAT  where  the  YY  is  the 
year  and  DDD  is  the  Julian  date.  The  .DAT  extension  signifies 
that  this  is  a  data  file.  For  this  thesis,  the  file  created 
was  90089MSIS.DAT.  Any  further  reference  to  this  file  will  be 
made  as  MSIS.  To  avoid  possible  confusion,  the  MSIS-83  model 
is  referred  to  as  MSIS3.  This  is  shortened  from  the  IDL 
version  of  MSIS-83  which  is  MSIS3.PRO. 

After  the  neutral  atmosphere  calculation  is  complete,  the 
user  is  asked  if  he/she  would  like  to  calculate  the 
photoel  ectron  ionization  rates  for  N2 ,  O2  and  0,  and  the 
dissociation  rate  for  N2 .  These  calculations  are  performed  by 
the  Naval  Research  Laboratory  (NRL)  PEGFAC  model  (Strickland 
and  Meier,  1982)  which  is  called  from  MSIS3  but  is  a  separate 
model.  This  creates  a  file  90089FLXM. PRT  where  the  M 
signifies  that  this  file  corresponds  to  the  MSIS  neutral 
atmosphere.   Further  references  to  90089FLXM. PRT  are  made  as 
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FLXM.   If  one  wishes  to  modify  the  neutral  atmosphere,  one 
must  use  the  next  model  called  MODATM. 

The  second  model,  MODATM,  allows  the  user  to  modify  the 
MSIS  neutral  atmosphere.  The  user  runs  MODATM  after  running 
MSIS.  MODATM  asks  for  the  date  in  IYD  format.  The  user  is 
able  to  change  the  following  inputs:  Tinf ,  T_120,  T_86,  O_120, 
O2_120,  N2_120  and  the  s  parameter.  The  s  parameter  is  used 
to  calculate  the  temperature  above  120  km.  This  parameter  is 
not  usually  changed.  Tinf  is  the  isothermal  temperature  at 
the  top  of  the  temperature  profile.  T_120  and  T_86  are  the 
temperatures  at  an  altitude  of  120  km  and  86  km  respectively. 
The  O_120,  O2_120,  and  N2_120  inputs  are  the  densities  for 
these  species  at  an  altitude  of  120  km.  When  modifying  the 
neutral  atmosphere,  the  MSIS  values  are  shown  as  a  reminder  to 
the  user.  For  this  thesis,  only  the  120  km  densities  and  Tinf 
were  changed.  Changing  T_86  and  T_120  has  much  less  affect  on 
the  profile  than  does  changing  Tinf.  MODATM  uses  the 
densities  at  120  km  to  re-calculate  the  density  profiles. 
This  re-calculation  is  affected  by  both  the  temperature  at  a 
given  altitude  and  Tinf.  At  lower  altitudes,  the  effects  of 
each  are  about  equal.  With  increasing  altitude,  Tinf 
dominates.  Therefore,  only  Tinf  was  adjusted  for  the  analysis 
presented  in  this  thesis.  Once  the  changes  have  been  made, 
MODATM  calculates  the  new  neutral  atmosphere.  It  stores  this 
information  in  a  file  called  YYDDDNEU"V" . DAT  where  the  "V"  is 
a  version  number  to  help  distinguish  between  various  runs. 
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For  example,  the  first  MODATM  neutral  atmosphere  file  for  this 
thesis  is  90089NEU0.DAT.  These  files  have  the  same  format  as 
the  MSIS  file. 

After  the  atmosphere  is  modified,  PEGFAC  is  again  used  to 
calculate  the  photoelectron  ionization  rates  and  the  N2 
dissociation  rate.  This  is  done  while  the  MODATM  is  being 
run.  This  information  is  stored  in  a  file  called 
YYDDDFLX"V".PRT  where  the  "V"  is  the  same  number  as  the  "V"  in 
the  corresponding  YYDDDNEUMV" .DAT  file.  Further  references  to 
the  MODATM  neutral  atmosphere  and  flux  files  will  be  made  as 
NEU"V"  and  FLX"V" .  Both  of  these  files  have  the  same  format 
as  the  MSIS3  files.  Now  that  the  calculation  of  a  neutral 
atmosphere  has  been  discussed,  the  procedure  for  calculating 
the  intensity  profiles  for  the  VK,  LBH,  and  OI2972A  emission 
must  be  explained. 

The  atmospheric  emission  model  is  a  collection  of  IDL 
procedures  that  calculates  the  intensity  profiles  for  the 
various  emissions.  The  IDL  codes  are  contained  in  the 
Appendix.  The  procedures  that  actually  calculate  the 
intensity  profiles  are  SLNTINT.PRO  and  OISLNT.PRO.  Although 
the  profiles  have  been  called  intensity  profiles,  they  are 
actually  slant  intensity  profiles.  A  slant  intensity  is  a 
column  intensity  at  a  given  viewing  angle.  Unless  specified, 
intensity  will  refer  to  a  slant  intensity. 

SLNTINT.PRO  calculates  the  intensities  for  VK  and  LBH. 
One  inputs  the  angle  between  the  zenith  and  the  instrument 
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viewing  direction.  This  accounts  for  the  difference  in 
pathlength  between  the  slant  and  vertical  column  intensities. 
For  a  detailed  description  of  this  see  (Cleary,  1985). 
SLNTINT  calls  a  procedure  called  GETVE.PRO.  GETVE.PRO  uses 
the  following  equation  to  calculate  the  volume  emission  rate 
for  the  (0,5)  VK  band: 


VVKs   (  (A3S+B3P+C3P/ 2  +  W3D)  xBO 0XQ05) 


*    ™  [O]  (2-16) 


In  equation  (2-16),  A3S,  B3P,  C3P,  and  W3D  are  the  excitation 
rates  for  the  A3E,  B3Tl,  CPU,  W3A  states,  respectively;  BOO  is 
the  branching  ratio  for  the  zeroth  vibrational  level  of  the  A 
state;  and  Q05  is  the  Frank-Condon  factor  for  the  (0,5) 
transition.  BOO  takes  into  account  that  the  A  state  is 
populated  by  direct  excitation  and  cascade  from  the  B,C,  and 
W  states.  The  second  term  in  the  denominator  accounts  for 
losses  in  the  A  state  due  to  quenching  by  atomic  oxygen.  QO 
is  the  oxygen  quenching  rate  and  A  is  the  Einstein  coefficient 
for  the  A  state.  Since  the  upper  state  for  LBH  is  populated 
entirely  by  excitation  to  the  a  state,  its  volume  emission 
rate  is  equal  to  the  A1P  rate.  The  rates  are  passed  into 
GETVE  from  the  desired  FLX"V"  file  by  the  procedure 
GETXRATE.PRO.  Densities  are  pass  into  GETVE  from  MSIS  or  the 
NEU"VM  file  by  GETDENS.PRO.  SLNTINT  then  uses  a  procedure 
called  VOLSLNT.PRO  to  convert  the  volume  emission  rates  to 
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slant  intensities.  SLNTINT  passes  out  an  altitude  array  and 
the  intensity  arrays  for  VK  and  LBH.  By  plotting  the 
intensity  arrays  versus  the  altitude  array,  the  intensity 
profiles  for  VK  and  LBH  are  produced.  Now  that  the  production 
of  the  VK  and  LBH  intensity  curves  has  been  discussed,  the 
production  of  the  curve  for  OI2972A  must  be  addressed. 

OISLNT.PRO  is  the  model  that  calculates  the  intensity  for 
OI2972A.  Prior  to  running  OISLNT,  one  must  run  NOX1DIM.PRO. 
This  is  a  one-dimensional  photochemical  model  which  calculates 
the  densities  for  C^*  and  N  S,  and  the  electron  density 
(Cleary,  1985).  The  0«'  and  electron  densities  are  stored  in 
a  file  called  YYDDDION"V" . DAT  and  the  N4S  is  stored  in  a  file 
called  YYDDDNOX"V" .DAT.  One  inputs  the  viewing  angle  and  a 
scaling  factor  into  OISLNT.  The  procedure  then  asks  the  user 
to  specify  the  desired  neutral  atmosphere  (i.e.,  MSIS,  NEUO , 
NEU1 ,  etc.).  Because  PCE  is  assumed,  OISLNT  calculates  the 
volume  emission  rate  by  dividing  the  production  terms  by  the 
loss  term  and  multiplying  this  by  A2072  •  The  first  production 
term  comes  from  the  reaction  given  by  equations  (2-5a)  and  (2- 
5b).  The  second  production  term  comes  from  equation  (2-7b). 
The  third  production  term  comes  from  equation  (2-7c).  The 
loss  term  is  the  sum  of  the  two  Einstein  coefficients.  (  To 
see  the  actual  equations  for  the  production  terms  see  the  code 
for  OISLNT.PRO  in  the  Appendix  ).  The  procedure  VOLSLNT 
converts  this  volume  emission  rate  to  a  slant  intensity.   The 
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resulting  intensity  array  is  plotted  versus  altitude  giving 
the  intensity  profile  for  the  OI2972A  emission. 

Now  that  the  background  information  has  been  covered,  the 
experiment  and  analysis  can  be  presented  starting  with  the 
experiment  in  the  next  chapter. 
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III.    THE  EXPERIMENT 

A.   MUSTANG  INSTRUMENT-CONSTRUCTION  AND  OPERATION 

The  MUSTANG  instrument  has  three  distinct  parts:  the 
telescope,  the  spectrograph,  and  the  electronics  interface. 
Figure  3-1  shows  a  drawing  of  the  MUSTANG  instrument.  A  brief 
description  of  each  section  and  their  function  follows. 

The  first  major  component  of  the  instrument  is  the 
telescope  section.  This  is  the  collecting  and  focusing  device 
for  the  instrument.  The  telescope  is  positioned  off-axis  from 
a  l/8th  meter  spherical  mirror.  Collimated  light  enters  the 
telescope  from  the  atmosphere  and  strikes  the  mirror  which 
focuses  the  light  at  the  entrance  of  the  spectrograph.  At 
this  point,  the  light  enters  the  next  section. 

The  second  section  of  the  instrument  is  a  modified  Ebert- 
Fastie  spectrograph.  A  5mm  x  140  urn  entrance  slit  at  the 
focus  of  the  telescope  mirror  acts  as  the  aperture  of  the 
spectrograph.  The  light  then  diverges  from  this  point  and 
falls  on  a  l/8th  meter  Ebert  mirror.  This  mirror  re- 
collimates  the  light  onto  the  reflective  diffraction  grating. 
This  1200  lines/mm  grating  causes  the  light  to  reflect  at  an 
angle  according  to  the  following  equation: 

asindm=mX.  (3-1) 
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Schematic  Diagram  of  MUSTANG  Instrument.  Major  components 
are:  (1)  1/8 th  m  Telescope  mirror,  (2)  Spectrograph  entrance  slit, 
(3)  l/8th  m  Ebert  mirror,  (4)  Diffraction  grating,  (5)  ITT  Image 
intensifier,  and  (6)  Hamamatsu  Image  detector. 


Figure   3-1   Schematic  drawing   of   the  MUSTANG   instrument 
(electronic  interface  not  shown) 
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The  light  then  reflects  off  the  Ebert  mirror  again  and  onto 
the  ITT  image  intensifier  in  the  1800-3400A  bandwidth.  Due  to 
the  diffraction  grating  and  the  mirror,  light  of  different 
wavelengths  falls  on  the  image  intensifier  in  different 
positions.  Once  calibrated,  it  is  the  position  that 
ultimately  allows  an  intensity/wavelength  relationship  to  be 
determined.  The  image  intensifier  is  comprised  of  a  quartz 
input  window,  cesium  telluride  (CsTe)  photocathode ,  two 
microchannel  plates  (MCP),  a  phosphor  coated  aluminum  screen, 
and  a  fiber  optic  output  window.  The  image  intensifier 
amplifies  and  converts  the  ultra-violet  light  to  visible 
light.  The  photocathode  generates  a  photoelectron  when  it  is 
struck  by  a  UV  photon.  The  electron  is  accelerated  by  an 
electric  potential  in  the  direction  of  the  microchannel 
plates.  When  the  electron  strikes  the  MCP ' s ,  an  electron 
avalanche  occurs.  The  MCP  acts  as  an  electron  multiplier 
producing  approximately  15,000  electrons  for  every  incoming 
electron.  Another  potential  difference  accelerates  these 
electrons  from  the  second  microchannel  plate  to  an  aluminum 
screen.  When  the  electrons  strike  the  phosphor  coating  they 
cause  the  screen  to  give  off  visible  light.  The  image 
intensifier  is  optically  coupled  to  a  Hamamatsu  linear  image 
sensor  by  its  fiber  optic  output  window  and  the  image  sensor's 
fiber  optic  input  window.  The  image  sensor  uses  512 
photodiodes  to  convert  the  visible  light  intensity  to  a 
proportional  voltage. 
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The  third  section  of  the  MUSTANG  is  the  electronic 
interface.  This  section  uses  an  analog-to-digital  converter 
and  a  first-in-first-out  (FIFO)  shift  register  to  send  the 
data  to  the  telemetry  section  of  the  rocket.  The  data 
provided  to  the  telemetry  is  sent  to  the  ground  via  a  RF 
signal.  For  a  detailed  description  of  the  interface  section 
see  (Quint,  1991). 

B.   THE  LAUNCH 

The  MUSTANG  instrument  was  launched  on  30  March  1990  at 
1000  local  time  (1700  GMT)  from  the  White  Sands  Missile  Range, 
New  Mexico.  The  launch  vehicle  was  a  NASA  Terrier  boosted 
Black  Brant  sounding  rocket.  Figure  3-2  shows  a  basic  diagram 
of  the  sounding  rocket.  The  MUSTANG  instrument  was  located  in 
the  experiment  section  of  the  payload.  Other  sections 
important  to  the  experiment  were  the  HV  section  which  provided 
the  power  source  and  the  T/M  section  which  transmitted  the 
data  back  to  the  ground.  The  payload  reached  an  apogee  of 
approximately  320  km.  Data  were  recorded  on  both  the  up-leg 
and  down-leg  of  the  flight.  This  experiment  resulted  in  the 
collection  of  airglow  data  in  the  form  of  8000  spectra  with  a 
wavelength  range  from  1800A  to  3400A.  The  spectra  were  taken 
between  an  altitude  of  100  and  320  km. 

In  order  to  determine  the  intensity  of  various  emission 
bands,  synthetic  spectra  of  these  emissions  were  fit  to  the 
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Figure  3-2  Launch  vehicle  block  diagram. 
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spectra  of  the  data.  Synthetic  spectra  are  spectra  calculated 
for  individual  emission  bands  by  models.  The  intensity  of 
these  spectra  can  be  adjusted.  By  adding  the  synthetic 
spectra  together  and  adjusting  their  intensities,  fits  for  the 
data  were  achieved.  Once  these  fits  were  obtained,  the 
intensities  for  the  various  emissions  were  assumed  to  be  the 
intensities  of  their  corresponding  synthetic  spectra. 

The  above  process  was  performed  by  Walden  (1991)  using  the 
N2  (0,5)  Vegard-Kaplan,  Nj  Lyman-Birge-Hopf ield,  2972A  atomic 
oxygen  and  several  other  emissions.  The  8000  spectra  were 
averaged  into  25  altitude  bins.  Each  bin  covered  a  10  km 
altitude  range.  As  a  result,  the  intensity  versus  altitude 
profiles  were  determined  for  these  emissions  for  the  1990 
flight.  The  profiles  for  the  VK,  LBH,  and  OI2972A  emissions 
were  used  in  the  atmospheric  emission  model  to  determine  the 
neutral  atmosphere. 

Now  that  the  background  information  and  the  experiment 
have  been  discussed  the  data  analysis  and  results  can  be 
addressed . 
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IV.     DATA  ANALYSIS  AND  RESULTS 

A.   ANALYSIS 

1.   Theory 

The  technique  used  to  analyze  the  data  is  simple.  One 
adjusts  the  inputs  to  the  model  until  a  satisfactory  fit  is 
obtained.  The  neutral  atmosphere  that  produces  this  fit  is 
assumed  to  represent  the  conditions  at  the  time  of  the 
observations.  In  this  thesis,  unless  otherwise  stated,  the 
data  refers  to  the  intensity  profiles  from  the  analysis 
performed  by  Walden  (1991). 

As  a  starting  point,  one  uses  the  MSIS3  model  to 
generate  a  neutral  atmosphere  and  the  atmospheric  emission 
model  to  calculate  an  intensity  profile  for  a  specific 
emission.  The  photoel ectron  excitation  rates  required  by  the 
model  to  produce  intensity  profiles  for  emission  bands  such  as 
the  Lyman-Birge-Hopf ield,  Vegard-Kapl an,  and  atomic  oxygen 
2972A  are  created  by  the  PEGFAC  and  stored  in  the  FLXM  data 
file.  One  then  compares  the  intensity  profile  of  the  data 
with  that  of  the  model.  If  the  two  curves  fit,  then  the 
atmospheric  densities  of  0,  Oj ,  and  Nj  (also  written  as  [O], 
[Oj],  and  [Nj])  are  the  same  as  those  calculated  by  MSIS3  and 
stored  in  the  MSIS  data  file. 
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If  the  model  calculation  does  not  agree  with  the  data, 
then  one  must  modify  the  neutral  atmosphere  by  adjusting  Tinf 
and/or  [0],  [02],  and  [Nj]  at  120  km  using  MODATM.  The  model 
is  again  used  to  generate  an  intensity  profile  and  the  effect 
that  these  changes  have  on  the  calculated  profile  is  noted. 
It  was  found  that  an  increase(decrease)  in  Tinf 
increases(decreases )  the  scale  height  of  the  calculated 
emission  while  increasing(decreasing)  [Nj]  shifts  the  curve 
up(down)  in  altitude  and  increasing(decreasing)  [0]  shifts  the 
profile  right(left).  It  was  also  found  that  changing  [Oj]  has 
little  effect  in  the  region  of  interest.  These  effects  are 
consistent  with  the  findings  of  Mack  (1991).  This  process  of 
adjusting  the  neutral  atmosphere  and  comparing  it  with  the 
data  could  lead  to  several  sets  of  inputs  which  fit  the  data 
reasonably  well  in  shape  but  not  necessarily  magnitude. 

Although  the  general  shape  of  the  profiles  may  agree, 
it  may  be  necessary  to  multiply  the  curve  by  some  constant  to 
fine  tune  the  fit.  This  is  called  a  scaling  factor.  Scaling 
factors  are  generally  smal 1 --usual ly  less  than  ten.  They  are 
significant  because  they  have  the  effect  of  increasing  the 
solar  EUV  flux  used  by  MSIS3  and  MODATM  to  predict  the  neutral 
atmosphere.  For  example,  a  scaling  factor  of  1.5  would 
increase  the  solar  EUV  flux  by  50%  which  would  have  the  end 
result  of  increasing  the  intensity  by  50%.  In  addition,  the 
same  scaling  factor  should  be  required  to  fit  all  other 
emission  bands  since  the  solar  flux  is  the  same  for  each.   An 
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inconsistency  in  the  scaling  factor  indicates  that  the  set  of 
inputs  used  to  produce  the  neutral  atmosphere  is  suspect.  The 
connection  between  intensity  and  the  solar  flux  was  described 
in  Chapter  1 1 . 

Once  a  set  of  likely  neutral  atmospheres  has  been 
determined,  then  another  emission  is  selected  for  examination. 
Intensity  profiles  for  this  emission  are  calculated  using  the 
atmospheric  emission  model  and  each  of  the  different  sets  of 
inputs.  The  profiles  for  this  second  emission  are  compared 
with  the  data  as  before.  Some  of  these  calculated  profiles 
will  fit  the  data  reasonably  well  while  others  will  have  the 
wrong  shape,  scale  height,  etc.  and  can  be  eliminated.  One 
continues  this  process  of  selecting  an  emission,  calculating 
the  intensity  curves  using  the  remaining  input  sets  and 
eliminating  ones  that  produce  curves  that  do  not  fit  the  data 
until  only  one  combination  remains.  The  density  profiles 
stored  in  the  NEU"V"  file  which  corresponds  to  these  initial 
conditions  are  those  of  the  neutral  species  for  that  day  and 
time.  The  density  profiles  of  the  various  species  of  ions, 
excited  species,  and  the  electrons  can  then  be  calculated 
using  NOX1DIM.PRO  and  the  corresponding  MSIS  or  NEU"V"  file. 
These  profiles  are  stored  in  90089ION"V" . DAT  ( 90089NOX"V" . DAT 
for  the  excited  species)  where  again  the  "V"  is  either  M  for 
MSIS  or  the  version  number  from  MODATM. 
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Now  that  a  general  description  of  the  theory  has  been 
discussed,  the  step-by-step  procedure  used  to  analyze  the  data 
can  be  presented. 
2 .   Procedure 

To  begin  the  analysis,  the  MSIS3  program  was  run  using 
the  following  values:  AP=69,  F10.7=186.5,  ave  F10.7=189.9  and 
x-ray  scaling  =  5  (Cleary,  private  communication,  1992).  This 
produced  a  Tinf=1255K  and  the  following  densities  at  120  km: 
[O]  =  8.8E10  cm"3,  [O2]  =  4.4E10  cm"3,  and  [N2]  =  3.2E11  cm"3.  It  also 
created  the  data  files  MSIS  and  FLXM.  Figure  4-1  is  a 
comparison  of  the  LBH  curve  generated  by  the  atmospheric 
emission  model  using  the  above  inputs  with  the  data.  As  one 
can  see,  the  calculated  intensity  profile  for  LBH  shows  a 
reasonable  fit  with  the  data  both  in  scale  height  and  in  shape 
above  approximately  120  km  when  scaled.  The  scaling  factor 
necessary  to  get  a  good  fit  is  1.8.  Although  the  model  gives 
profile  information  below  120  km,  no  attempt  to  fit  the 
prediction  to  the  data  was  made  at  lower  altitudes  for  the 
following  reason.  Below  120  km,  it  is  difficult  to 
distinguish  between  the  LBH  bands  and  the  brighter  NO 
emissions.  This  caused  the  fit  between  the  spectra  and  the 
data  to  be  less  certain.  The  intensity  profile  for  the  data 
obtained  from  the  fit  in  the  region  below  120  km  is  therefore 
less  reliable. 
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Figure  4-1   Comparison  of  the  data  with  the  model  intensity 
profile   for  LBH  using  MSIS  values:   [O]_120=8 . 8E10   cm" 
[O2]_120  =  4.4E10  cm"3,  [N2]_120  =  3 .  2E11  cm"3  and  Tinf  =  1255K. 
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Even  though  the  above  trial  produced  a  good  fit,  two 
other  combinations  of  inputs  were  tried.  These  sets  were 
suggested  by  Mack  (1991).  The  purpose  of  this  was  to  test 
whether  either  of  these  two  neutral  atmospheres  would  lead  to 
a  fit  when  compared  with  the  data  from  Walden's  analysis.  To 
do  this,  MODATM  was  used  to  change  the  value  of  Tinf  to  1500K 
and  to  reduce  the  MSIS3  [0]  at  120  km  by  27%  and  increase  the 
MSIS3  [02]  and  [N2]  at  120  km  by  5%  and  30%  respectively.  This 
trial  produced  NEUO  and  FLXO .  These  changes  were  made  based 
on  values  suggested  on  page  82  of  the  report  by  Mack  (1991). 
The  effects  of  these  modifications  can  be  seen  by  the 
comparison  of  the  LBH  intensity  curves  of  the  atmospheric 
emission  model  with  the  data  shown  in  Figure  4-2.  The  curve 
produced  from  these  values  has  the  wrong  scale  height  and  it 
also  required  scaling  which  was  not  needed  by  Mack. 

The  other  possible  combination  used  Tinf=1700K  and  the 
input  values  for  the  neutral  species  densities  at  120  km  as 
suggested  on  page  83  of  Mack  (1991).  These  values  represent 
a  48%  reduction  in  [O],  a  6%  increase  in  [O*]  and  the  MSIS3 
value  for  [Nj]  at  120  km.  A  run  using  these  values  was 
performed  and  the  calculations  were  stored  in  NEU1  and  FLX1 . 
Figure  4-3  shows  the  intensity  profile  generated  by  the 
atmospheric  emission  model  for  LBH  produced  by  this 
modification  and  the  data  profile.  One  can  see  that  ^the  curve 
produced  using  these  inputs  also  does  not  agree  wi1fh~  the  data. 
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Figure  4-2  Comparison  of  the  data 
profile  for  LBH  using  MODATM  values: 
+  5%  [O2]_120  and  +30%  [N,  ]_120. 


with  the  model  intensity 
Tinf=1500K,  -27%  [O]_120, 
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Figure  4-3  Comparison  of  the  data  with  the  model  intensity 
profile  for  LBH  using  MODATM  values:  Tinf=1700K,  -48%  [0]  120 
and  +6%  [Oj]_120.  ~ 
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When  compared  with  the  intensity  profiles  for  the  data  from 
(Walden,  1991),  neither  set  of  inputs  fit  the  data.  The  fact 
that  this  analysis  rejected  the  sets  of  possible  inputs 
suggested  by  Mack  (1991)  can  be  explained  in  the  following 
way.  The  intensity  profiles  for  the  data  were  determined  by 
fitting  synthetic  spectra  to  the  data  spectra.  Mack  (1991) 
used  only  part  of  the  MUSTANG  wavelength  range  while  the  later 
analysis  by  Walden  (1991)  used  the  full  wavelength  range  of 
the  instrument.  Analysis  over  a  smaller  wavelength  range, 
which  includes  fewer  emission  features,  is  much  more  sensitive 
to  large  fluctuations  in  the  intensity  than  analysis  over  a 
larger  range.  One  can  imagine  that  if  a  fit  was  performed  at 
the  only  large  data  spike  in  an  otherwise  even  spectra,  the 
fit  would  be  poor  for  the  whole  spectral  range.  Because  of 
this,  the  data  profiles  produced  by  Walden  (1991)  are 
considered  more  accurate.  The  input  sets  suggested  by  Mack 
(1991)  do  not  fit  the  improved  data  curves.  Therefore  they 
were  rejected. 

Since  the  original  curve  from  MSIS/FLXM  showed  better 
agreement  than  either  of  the  profiles  from  NEUO/FLXO  and 
NEU1/FLX1,  an  attempt  was  made  to  fit  the  data  using  the 
original  Tinf=1255K  while  varying  the  neutral  densities. 
Pursuit  of  a  fit  using  NEUO/FLXO  or  NEU1/FLX1  was  not 
continued  because  the  scale  heights  of  these  curves  did  not 
match  the  data.   It  was  apparent  that  their  Tinf  values  were 
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not  correct  and  that  the  modification  of  the  densities  would 
not  result  in  a  fit. 

Continuing  the  analysis,  MODATM  and  the  original  MSIS3 
values  for  Tinf,  [O2]_120,  and  [N2]_120  were  used.  Only  half 
of  the  [O]_120  was  used  and  NEU2/FLX2  was  created.  From 
Figure  4-4  which  shows  this  LBH  intensity  profile  and  the 
data,  one  can  see  that  this  modification  produced  a  curve  with 
a  peak  intensity  which  was  a  little  too  low.  As  stated 
earlier,  increasing  [N2]_120  causes  the  peak  intensity  of  the 
model  calculation  to  occur  at  a  higher  altitude.  The  next 
trial  makes  use  of  this  effect.  To  create  the  files  NEU3  and 
FLX3,  MODATM  was  run  with  all  NEU2/FLX2  values  except  [N2]_120 
was  increased  by  30%.  The  intensity  profile  produced  with 
these  values  is  shown  in  Figure  4-5.  The  effects  of  this 
change  can  be  seen  by  comparing  Figures  4-4  and  4-5.  This 
last  run,  when  scaled  by  1.3,  showed  excellent  agreement 
between  the  model  and  the  data.  At  this  point,  two  input 
sets,  MSIS  and  NEU3 ,  still  look  capable  of  producing  the 
correct  neutral  atmosphere  so  it  was  necessary  to  choose 
another  emission  and  to  continue  the  analysis. 

The  next  emission  profile  examined  was  the  (0,5) 
Vegard-Kaplan  band.  When  the  model  intensity  profiles  using 
MSIS  and  NEU3  were  compared  with  the  data,  neither  of  the 
profiles  for  VK  agreed  well  with  the  data  as  seen  in  Figures 
4-6  and  4-7.  Both  curves  fit  the  peak  intensity  of  the  data 
well.  However,  moving  above  the  peak,  the  agreement  for  both 
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Figure  4-4 
prof i le  for 


Comparison  of  the  data  with  the  model  intensity 
LBH  using  MODATM  values:  Tinf  =  1255K  and  -50%  [0]  . 
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Figure  4-5  Comparison  of  the  data  with  the  model  intensity 
profile  for  LBH  using  MODATM  values:  Tinf=1255K,  -50%  [O]_120, 
and  +30%  [N,]_120. 
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Figure  4-6  Comparison  of 
profile  for  VK  using  MSJS 
[O2]_120  =  4.4E10  cm"', 


cm 


the  data  with  the  model  intensity 
values:  Tinf=1255K,  [O]_120=8 . 8E10 
and  [N,]_120  =  3.2E11  cm*5. 
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Figure  4-7   Comparison  of  the  data  with 
profile  for  VK  using  MODATM  values: 
[N2]_120  and  Tinf=1255K. 


the  model  intensity 
-50%   [O]_120,   +30% 
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worsened  due  to  a  difference  in  the  scale  heights. 
Furthermore,  in  order  to  get  the  fit  at  the  peak,  it  was 
necessary  to  scale  the  curve  for  MSIS  by  4.8  and  2.7  for  NEU3. 
Because  the  solar  flux  is  the  same  for  all  emissions,  the 
scaling  factor  should  be  the  same  regardless  of  the  emission 
being  examined.  Therefore,  any  large  variation  in  the  scaling 
factor  is  an  indicator  that  that  set  of  inputs  should  be 
considered  suspect. 

Because  neither  MSIS  nor  NEU3  lead  to  a  good  fit,  the 
order  in  which  the  emission  features  were  fitted  was  changed. 
The  intensity  profile  of  the  VK  band  is  more  reliable  data 
than  the  profile  for  the  LBH  because  the  signal-to-noise  ratio 
of  the  VK  data  was  greater.  As  a  result,  reversing  the  order 
of  the  analysis  to  VK  and  then  LBH  is  a  better  approach. 
Furthermore,  this  ordering  is  superior  because  the  model  for 
VK  is  more  complicated  due  to  a  greater  number  of  production 
terms.  This  causes  more  variables  to  be  constrained  by 
analyzing  the  VK  bands  first. 

To  begin  the  analysis  using  VK,  a  value  of  3E-11  for 
the  atomic  oxygen  quenching  rate  coefficient  was  determined  by 
taking  the  average  of  the  values  reported  by  Piper  et  al  . 
(1981),  Thomas  and  Kaufman  (1985)  and  De  Souza  et  al .  (1985). 
Also,  an  estimate  of  0.125  for  the  zero  vibrational  level 
branching  ratio  was  used.  This  was  determined  in  the 
following  manner.  The  first  8  vibrational  levels  (0  through 
7)  for  A3Z  (no  quench)  from  Figure  22  of  (Cartwright,  1978) 
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appear  to  be  of  the  same  order  of  magnitude.  Therefore,  the 
branching  ratio  for  the  zeroth  vibrational  level  was  estimated 
to  be  1/8  of  the  total  relative  number  density  of  these  states 
or  0.125.  Using  these  values  in  the  atmospheric  emission 
model,  previously  rejected  runs  were  re-evaluated.  NEU0/FLX0 
and  NEU2/FLX2  were  again  rejected  because  the  curves  showed 
poor  agreement  with  the  data  (see  Figures  4-8  and  4-9).  On 
the  other  hand,  the  run  using  NEU1/FLX1  fit  the  data  well 
except  for  a  slight  difference  in  the  scale  height  which  was 
only  a  little  too  high  as  seen  in  Figure  4-10.  Although  the 
calculations  using  MSIS/FLXM  and  NEU3/FLX3  were  found  to  be 
poor  as  stated  above,  they  were  the  best  two  fits  from  the 
original  analysis  attempt.  Therefore,  they  were  not 
immediately  rejected.  In  addition,  these  runs,  like  the 
NEU1/FLX1  trials,  fit  pretty  well  at  the  peak  intensity  but 
differed  from  the  data  elsewhere  due  to  the  scale  height 
difference . 

One  way  in  which  the  scale  height  can  be  increased  is 
to  increase  the  atomic  oxygen  quenching  rate  coefficient. 
Figure  4-11  shows  a  VK  intensity  curve  using  3E-11  and  one 
using  9E-11.  This  illustrates  the  effect  that  this  change  has 
on  the  VK  intensity  profile.  This  higher  coefficient  was 
suggested  by  Sharp  (1971).  Unfortunately,  a  change  in  the 
scale  height  is  not  the  only  effect  that  increasing  the  atomic 
oxygen  quenching  rate  has  on  the  VK  profile.  Figure  4-11  also 
shows  that  by  increasing  the  rate,  the  peak  intensity  occurs 
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Figure  4-8   Comparison  of  the 
profile  for  VK  using  MODATM  va 


+  5%  [0-]_120  and  +30%  [N-  ]_120. 


data  with  the  model  intensity 
ues:  Tinf=1500K,  -27%  [O]_120, 
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Figure  4-9  Comparison  of  the  data  with  the  model  intensity 
profile  for  VK  using  MODATM  values:  Tinf=1255K  and  -50% 
[O]_120. 
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Figure  4-10  Comparison  of  the  data  with  the  model  intensity 
profile  for  VK  using  MODATM  values:  Tinf=1700K,  -48%  [O]_120 
and  +6%  [N2]_120. 
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Figure  4-11   Effect  of  changing  the  atomic  oxygen  quenching 
rate  coefficient  on  VK  intensity  profile. 
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at  a  higher  altitude  and  the  overall  intensity  drops.   This 
drop  in  intensity  is  most  pronounced  at  the  peak.   The 
MSIS/FLXM,  NEU1/FLX1,  and  NEU3/FLX3  runs  were  repeated  using 
the  higher  quenching  rate  coefficient  but  the  fits  were  again 
poor  and  they  required  an  even  larger  scaling  factor. 

Another  way  to  increase  the  scale  height  of  the  VK 
emission  profile  is  to  increase  Tinf  as  was  discussed  earlier. 
Taking  this  into  consideration,  one  might  guess  that  a 
slightly  lower  Tinf  might  improve  the  fit  using  NEU1/FLX1. 
Several  attempts  were  made  using  Tinf=1600K,  a  quenching  rate 
coefficient  of  9E-11  and  various  combinations  of  neutral 
densities  at  120  km.  These  trials  were  all  variations  of 
NEU4/FLX4  and  none  produced  the  desired  effects  (see  Figures 
4-12  through  4-14).  Scaling  was  attempted  to  fit  the  peak. 
Next,  the  quenching  rate  coefficient  was  changed  back  to  3E-11 
and  better  results  were  achieved  as  one  can  see  from  Figure  4- 
15.  In  this  run,  the  peak  intensity  was  too  high  so  NEU5/FLX5 
was  created  using  Tinf=1600k,  half  the  MSIS3  [O]_120  value  and 
increasing  the  MSIS3  [N2J_120  by  30%.  This  produced  better 
scale  height  agreement  but  the  altitude  for  the  peak  intensity 
was  still  a  little  too  high  as  one  can  see  from  the  intensity 
profile  in  Figure  4-16.  The  final  run  for  VK  used  the 
NEU5/FLX5  values  except  the  peak  was  lowered  by  lowering  the 
[N2]_120  back  to  the  original  MSIS3  [N2]_120.  This  run  was 
called  NEU6/FLX6  and  gave  an  excellent  fit  when  scaled  by  3.25 
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Figure  4-12  Comparison  of  the  data  with  the  model  intensity 
profile  for  VK  using  MSIS  values  and  changing  the  atomic 
oxygen  quenching  rate  coefficient  to  q^=9E-ll. 
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Figure  4-13   Comparison  of  the  data  with  the  model  intensity 
le  for  VK  using  MODATM  values:  Tinf=1600K,  -50%  [O]_120 

anging  to  q*=9E-ll. 
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Figure  4-14  Comparison  of  the  data  with  the  model  intensity 
profile  for  VK  using  MODATM  value  Tinf=1600  and  changing  to 
qft=9E-ll. 
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Figure  4-15 
profile  for 
q,=3E-ll. 


Comparison  of  the  data  with  the  model  intensity 
VK  using  MODATM  value  Tinf=1600K  and  changing  to 
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Figure  4-16  Comparison  of  the  data  with  the  model  intensity 
profile  for  VK  using  MODATM  values:  Tinf=1600K,  -50%  [O]_120, 
+30%  [N,]_120  and  changing  to  q,  =  3E-ll. 
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(see  Figure  4-17).  This  set  of  parameters  was  then  used  to 
fit  the  LBH  data.  This  trial  produced  a  good  fit  for  the  LBH 
only  at  the  peak  intensity  and  required  scaling  the  prediction 
by  1.5  (see  Figure  4-18).  A  good  fit  over  the  full  range  of 
the  profile  could  not  be  achieved.  Due  to  the  nature  of  the 
log/linear  plot,  any  difference  between  the  curve  and  the  data 
will  be  greatest  at  the  peak.  Therefore,  the  best  fit  was 
determined  to  be  the  one  that  best  matched  the  peak. 
Additionally,  the  signal  is  the  largest  at  the  peak  so  the 
spectral  analysis  that  produced  the  data  points  is  most 
reliable  there.  The  difference  between  the  scaling  factor 
required  for  the  VK  fit  and  the  LBH  fit  is  addressed  below. 

Previously,  the  analysis  of  VK  was  performed  using  an 
estimate  for  the  branching  ratio  ,  denoted  BOO.  A  more 
accurate  value  for  BOO  was  determined  by  taking  the  ratio  of 
the  relative  number  density  at  the  zero  vibrational  quantum 
number  to  the  total  relative  number  density  for  the  N2  excited 
state  A.  The  new  value  for  BOO  was  found  to  be  0.25  and 
Figure  4-19  shows  the  effect  that  this  change  has  on  the  VK 
intensity  profile.  This  curve  required  a  scaling  factor  of 
1.6  to  achieve  a  good  fit.  This  is  consistent  with  the 
scaling  factor  of  1.5  required  for  LBH  curve  as  shown  in 
Figure  4-18. 

The  final  emission  used  to  analyze  the  data  was  the 
OI2972A  band.   Using  the  model  procedure  OISLNT.PRO  and 
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Figure  4-17  Comparison  of  the  data  with  the  model  intensity 
profile  for  VK  using  MODATM  values:  Tinf=1600K,  -50%  [O]_120 
and  changing  to  q,»=3E-ll. 


57 


I              ■   ■    T 1 

1 

1 

1— 

i 

& 

- 

**\ 

pi 

X 

- 

< 

1— 

. 

. 

< 

* 

X 

Q 

1 

- 

X 

- 

1— 

^ 

/    * 
/   * 

t 

t/i 

/    * 

* 

z 

LU 

/    * 

K 

Z 

/         X 

_ 

X 

1     * 

CD 

_J 

/     * 

U_ 

• 

O 

in 

- 

Z 

/      * 

*-' 

o 

en 

/       * 

>- 

CD 

/v 

O 

■ 

< 

UJ 

_J 

Q_ 

/        * 

< 

^ 

CO 

o 
o 

>- 
or 

* 

o 

UJ 

X 

• 

* 

* 

* 

1              1              1 

1 

1 

' 

i 

If) 
o 


o 
o 


o 


CO 

I 
o 

UJ 

_J 

cc 


10 

z 

UJ 

I- 


o 
o 

CM 


o 

O       T- 


(W>0  3QfUll1V 


Figure  4-18  Comparison  of  the  data  with  the  model  intensity 
profile  for  LBH  using  MODATM  values:  Tinf=1600K  and  -50% 
[O]_120. 


58 


■     1                       ! 

1 

i 

i 

_ 

- 

- 

- 

K        * 

- 

- 

J** 

^* 

- 

- 

r    * 

- 

- 

(0 

- 

>- 

m 

/ 

Q 
UJ 
_J 
< 
C_> 
if) 

- 

/ 

O 

UJ 

X 

- 

/                 1                       1 

i      i      i 

1 

* 

♦ 
* 

• 

o 
o 


o 
o 
o 

-I   o 


o 
o 
o 


o 

UJ 

_l 
>• 
< 


if) 

UJ 


o 
o 

CM 


O 

o 


(m)  aaruinv 


Figure  4-19   Comparison  of  the 
profile  for  VK  using  MODATM 
[O]_120;  q,=3E-ll  and  changing 


data  with  the  model  intensity 
values:   Tinf=1600K  and  -50% 
to  B00=0.25. 
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NEU6/FLX6,  a  predicted  intensity  profile  was  produced  for  this 
emission  and  compared  with  the  data.  As  discussed  in  Chapter 
II,  the  profile  calculated  by  OISLNT.PRO  is  only  accurate  for 
altitudes  above  115  km.  Because  of  this,  only  the  upper 
altitudes  were  considered  when  attempting  to  fit  the  model 
profile  to  the  data.  A  fair  fit  was  achieved  when  the  curve 
was  scaled  by  2.1.  This  is  shown  by  the  comparison  of  the 
model  OI2972A  intensity  profile  with  the  data  in  Figure  4-20. 
Even  though  the  analysis  for  the  OI2972A  emission  produced 
only  a  fair  fit,  it  is  important  to  remember  that  this  is  the 
least  reliable  band.  Because  a  problem  with  the  telemetry 
caused  some  data  to  be  lost  in  the  region  of  this  wavelength 
(Walden,  1991),  there  is  a  systematic  uncertainty  in  the  2972A 
intensity  profile.  With  the  systematic  error  and  the  error 
inherent  in  the  manual  fitting  technique,  the  2.1  scaling 
factor  is  acceptable. 

The  above  analysis  indicated  that  the  parameters  from 
NEU6/FLX6  should  correspond  to  the  neutral  atmosphere  at  the 
time  of  the  launch.  To  test  this,  NOX1DIM.PRO  (Cleary,  1986) 
was  then  run  using  these  files  and  the  profile  for  the 
electron  density,  [ e" ] ,  was  produced.  This  procedure  uses  the 
one  dimensional  MCE.  A  discussion  of  the  results  and  a 
comparison  of  the  calculated  [e"]  profile  with  that  of  an 
ionosonde  sounding  performed  by  the  Geophysics  Laboratory  (GL) 
at  the  time  of  the  launch  are  presented  in  the  following 
section.   No  comparison  was  made  for  the  charged  species 
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Figure  4-20  Comparison  of  the  data  with  the  model  intensity 
profile  for  OI2972A  using  MODATM  values:  Tinf=1600K  and  -50% 
[O]_120. 
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because  no  independent  density  profiles  of  these  species  exist 
for  this  launch. 

In  the  course  of  writing  this  thesis,  it  was  discovered 
that  the  W3D  production  rate  was  inadvertently  omitted  from 
the  calculation  for  VK  (see  Equation  2-16).  This  omission 
caused  a  larger  scaling  factor  to  be  required  for  all  fits  of 
VK .  When  the  term  was  added,  it  was  found  that  the  best  fit, 
NEU6/FLX6,  fit  without  scaling.  This  is  shown  in  Figure  4-21. 
Figures  4-22  and  4-23  shows  two  other  figures  for  comparison 
to  show  that  the  ratio  of  the  old  scaling  factor  to  the  new 
scaling  factor  is  approximately  1.6.  This  shows  that  the 
missing  term  was  incorporated  in  the  scaling  factor.  This 
applies  only  to  VK  and  an  approximate  1.5  scaling  factor  is 
still  required  for  LBH  and  012972.  This  violates  the 
requirement  that  the  scaling  factor  for  each  emission  be  the 
same.  If,  for  example,  the  scaling  was  accounting  for  a  value 
for  the  solar  flux  which  was  too  low,  then  the  scaling  factor 
would  have  to  be  the  same  for  all  emissions.  The  scaling 
factor  is  still  important  as  an  indicator  of  goodness  of  fit. 
Because  the  omission  of  the  W3D  term  only  resulted  in  a 
magnitude  difference  in  all  the  VK  plots,  the  above  analysis 
would  still  have  led  to  NEU6/FLX6  as  the  best  fit. 
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Figure  4-21 
production. 


Re-plot   of  Figure  4-19  using  correct   VK 
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Figure  4-22 
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Re-plot   of  Figure   4-6  using  correct  VK 
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Figure  4-23 
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Re-plot   of  Figure  4-15  using  correct   VK 
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B.   RESULTS 

The  above  analysis  lead  to  an  estimate  of  Tinf=1600K  and 
a  value  for  atomic  oxygen  density  at  120  km  of  4.4E10  cm 
which  is  half  the  MSIS3  calculated  value.  If  NEU6  does  in 
fact  correspond  to  the  neutral  atmosphere  at  the  time  of  the 
launch,  then  the  [ e" ]  curve  calculated  by  NOX1DIM  using  NEU6 
would  be  expected  to  fit  the  [e"]  profile  measured  by  GL . 
Figure  4-24  compares  the  NOX1DIM.PRO  calculated  [e*]  profile 
using  NEU6  with  the  [e~]  profile  from  the  GL  ionosonde 
measurement  (Bullett  and  Buchau,  1990).  Figure  4-25  compares 
two  NOX1DIM.PRO  [ e"  ]  calculations  with  the  ionosonde  data:  one 
using  the  MSIS  neutral  atmosphere  and  the  other  using  the  NEU6 
neutral  atmosphere.  This  is  intended  to  show  whether  any 
improvement  can  be  seen  in  the  calculated  electron  density 
profile  by  using  the  method  describe  above  to  determine  the 
neutral  atmosphere.  As  one  can  see,  neither  calculation 
agrees  well  with  the  GL  data.  The  MSIS  and  the  MODATM  curves 
do  seem  to  have  the  same  general  shape  as  the  data  from  90  to 
160  km  and  90  to  280  km  respectively.  The  conclusions  drawn 
from  the  above  results  follow  in  the  next  chapter. 
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Figure  4-24  Comparison  of  the  ionosonde  data  with  the  N0X1DIM 
calculated  [e"]  profile  using  MODATM  values:  Tinf=1600K  and 
-50%  [O]_120. 
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Figure  4-25  Comparison  of  the  ionosonde  data  with  the  N0X1DIM 
calculated  [e']  profiles  using  MSIS  values  and  using  NEU6/FLX6 
val ues . 


68 


V.    CONCLUSIONS 
A .   SUMMARY 

It  was  the  aim  of  this  thesis  to  develop  a  new  technique 
for  determining  the  densities  of  0,  O2  and  N2  in  the 
ionosphere.  These  species  make-up  the  neutral  atmosphere.  It 
was  hoped  that  by  comparing  calculated  intensity  profiles  with 
the  data,  one  could  adjust  the  inputs  until  the  model 
intensity  profiles  fit  the  data.  The  neutral  atmosphere 
calculated  by  the  model  would  then  correspond  to  the  neutral 
atmosphere  at  the  time  that  the  data  was  taken.  This  could 
then  be  used  to  determine  the  densities  of  the  charged  species 
and  electrons. 

The  analysis  method  presented  in  this  thesis  did  not,  at 
first  glance,  appear  to  provide  a  viable  approach  for 
determining  the  neutral  atmosphere  or  the  densities  of  the 
charged  species  and  electrons.  This  is  somewhat  misleading 
for  several  reasons.  First,  the  data  from  the  GL  ionosonde 
measurement  was  taken  310  km  north  of  the  White  Sands  launch 
site.  Normally,  this  would  provide  an  accurate  electron 
density  profile  for  White  Sands  but  at  the  time  of  the  launch 
there  was  higher  than  normal  geomagnetic  activity.  This 
resulted  in  a  Traveling  Ionospheric  Disturbance  or  TID  which 
distorted  the  ionosphere  causing  lower  peak  densities  to  occur 
at  higher  altitudes.   These  effects  were  reported  by  Bullett 
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and  Buchau  (1990)  but  could  not  be  taken  into  account  in  the 
models  used  in  the  analysis.  Also,  since  the  ionosonde 
measurement  was  not  taken  in  White  Sands,  the  effects  of  the 
TID  at  the  launch  site  could  only  be  estimated  and  thus  the 
[e"]  is  also  an  estimate.  Second,  the  experiment  lost  some 
data  in  the  2972A  wavelength  region  due  to  a  problem  with  the 
telemetry  circuit  (Walden,  1991).  Third,  only  the  three 
emission  bands  were  used.  Using  other  prominent  emissions 
would  have  helped  the  analysis.  One  emission  that  would  have 
been  particularly  useful  is  the  Second  Positive  band  of 
molecular  nitrogen.  An  analysis  of  this  prominent  band  could 
not  be  performed  because  a  synthetic  spectrum  model  for  the 
data  currently  does  not  exist.  Another  method  for  determining 
the  intensity  profile  for  the  Second  Positive  suggested  by 
Walden  (1991)  proved  to  be  unworkable.  Had  the  OI2972A 
emission  been  more  reliable  and  if  the  Second  Positive  band 
had  been  available,  a  better  calculation  might  have  been 
possible  using  the  above  analysis  method.  Finally,  the 
scaling  factor  of  approximately  1.6  which  was  necessary  to 
achieve  a  good  fit  for  LBH  and  OI2972A  is  not  understood.  It 
suggests  that  there  is  uncertainty  in  the  models  for  the  LBH 
and  OI2972A  emissions.  The  excellent  fit  of  VK  also  suggests 
that  the  data  for  both  LBH  and  OI2972A  is  less  reliable  than 
originally  believed.  One  last  conclusion  that  can  be  drawn 
from  this  thesis  is  that  the  excellent  fit  for  VK  and  LBH 
(once  scaled)  in  scale  height  indicates  a  Tinf=1600K  which  is 
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higher  than  the  Tinf  calculated  by  the  MSIS3  model.  It  also 
is  different  than  the  1500K  and  1700K  suggested  by  Mack 
(1991).  The  Tinf=1600K  resulted  in  much  better  scale  height 
agreement  for  the  intensity  curves. 

The  method  described  in  this  thesis  did  not  produce  a 
neutral  atmosphere  which  led  to  an  electron  density  profile 
that  fit  the  Geophysics  Laboratory  data.  Although  this  is 
true,  this  approach  cannot  be  dismissed.  Improving  the 
problem  areas  discussed  above  should  allow  this  method  to 
successfully  calculate  the  neutral  atmosphere  as  well  as  the 
charged  species  and  electron  density  curves. 

B.   SUGGESTIONS  FOR  FOLLOW-ON  RESEARCH 

As  pointed  out  in  the  summary  section,  there  are  several 
areas  in  which  improvements  can  be  made  to  better  test  the 
viability  of  this  method  for  analyzing  the  ionosphere.  First, 
improvements  were  made  for  the  1992  launch  which  should  have 
increased  the  accuracy  of  the  data  collected.  One  such 
improvement  was  better  telemetry  electronics  which  eliminated 
the  data  drop  problems  that  hurt  the  information  collection  in 
1990.  Using  this  better  data,  the  synthetic  spectrum  fitting 
performed  by  Walden  (1991)  should  be  repeated  to  get  the 
intensity  profiles  for  VK,  LBH  and  OI2972A.  The  elimination 
of  the  drops  in  the  data  should  result  in  a  more  accurate 
intensity  profile  for  OI2972A.  Additionally,  the  model  for 
OI2972A  should  be  extended  to  include  the  lower  altitude 
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profile  of  this  emission.  Second,  a  synthetic  spectrum  should 
be  developed  for  the  Second  Positive  band.  As  was  suggested 
by  Walden  (1991),  the  known  contributions  of  Nj  could  be 
subtracted  from  the  data  and  the  remainder  could  be  assumed  to 
be  the  Second  Positive  spectrum.  This  should  now  be  practical 
with  the  better  1992  launch  data.  Third,  other  emissions 
should  also  be  used  in  the  fit  process.  Using  several 
emissions  should  lead  to  better  results.  Fourth,  the  scaling 
requirement  should  be  more  fully  investigated.  Possible 
refinements  in  the  models  for  LBH  and  012972  as  well  as  the 
better  data  for  the  1992  launch  should  be  able  to  eliminate 
these.  Lastly,  an  effort  should  be  made  to  fully  integrate 
the  atmospheric  emission  model,  the  MSIS3  model,  MODATM,  and 
N0X1DIM  into  a  single  program  that  uses  numerical  methods  to 
automatically  determine  the  fit.  This  would  be  necessary  to 
make  this  method  practical  for  analyzing  the  large  amount  of 
data  received  from  a  satellite-based  spectrograph.  If  the 
above  changes  are  implemented,  this  method  and  the  data  from 
satellites  could  be  used  to  make  timely  ionospheric  profile 
calculations  possible  for  any  point  on  the  globe. 
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APPENDIX 

The  Appendix  contains  the  IDL  code  for  the  Atmospheric 
Emission  Model . 


;THIS  PROCEDURE  GETS  THE  0,02,N2 

/DENSITIES 

;AND  CORRESPONDING  ALT  FROM  EITHER 

;MSIS.DAT  OR  NED'V'.DAT 

; FILES  FOR  USE  IN  GETVE.PRO. 

/WRITTEN  BY  KENNETH  L.  MCELROY 

;15  MAY  1992 

PRO  GETDENS , OALT, ODEN, 02DEN, N2DEN, TEMP 

FILE-'  ' 

READ, 'WHAT  DATA  FILE  DO  YOU  WANT  TO  USE?  (I.E.  9008  9MSIS.DAT  OR  9008 9NEU"V" . DAT) ' 

,FILE 

OPENR,2,FILE 

B-ASSOC(2,FLTARR(183) ) 

OALT-B(l) 

TEMP-B  (2) 

ODEN-B(3) 

D2DEN-B(4) 

N2DEN-B(5) 

CLOSE, 2 

RETURN 

END 
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;THIS  PROG  READS  THE  ALT  AND  VOL 
;PROD  RATES  A1P, A3S, B3P, C3P,  AND  W3D 
;FOR  N2  FROM  FLX'V'.PRT  OR 
;FLXM.PRT  FILES 
; WRITTEN  BY  K.L.  MCELROY 
;23  APR  1992 
/MODIFIED  06  JUNE  92 
PRO  GETXRATE,ALT,A1P, A3S,B3P,C3P, W3D,OI 

file-'  ' 

read, 'What  file  do  you  want  to  use? (FLX"V" .PRT  OR  FLXM . PRT) ', file 

OPENR, 1,FILE 

;THIS  SECTION  GETS  ALT  AND  A3S 
;FOR  N2 
STR1-'  ' 

WHILE  STRPOS (STR1, 'VOLUME  PRODUCTION  RATES')  EQ  -1  DO  READF,1,STR1 
FOR  J-1,5  DO  READF,1,STR1 
ALT-FLTARR(51) 
A3S-FLTARR(51) 
Rl-FLTARR(ll) 
FOR  J-0, 50  DO  BEGIN 

READF, 1,R1 

ALT (J) -R1  (0) 

A3S  (J)-Rl  (10) 
ENDFOR 

;THIS  SECTION  GETS  A1P,B3P,C3P 
;AND  W3D 
STR2-'  ' 

FOR  J-1,7  DO  READF,  1,STR2 
A1P-FLTARR(51) 
B3P-FLTARR(51) 
C3P-FLTARR(51) 
W3D-FLTARR(51) 
R2«FLTARR(5) 
FOR  J-0, 50  DO  BEGIN 

READF, 1,R2 

A1P (J)-R2 (4) 

B3P (J)-R2 (1) 

C3P (J)-R2 (2) 

W3D(J)-R2  (3) 
ENDFOR 
WHILE  STRPOS (STR2, 'OI' )  EQ  -1  DO  READF, 1,STR2 
FOR  J-1,3  DO  READF, 1, STR2 
OI-FLTARR(51) 
FOR  J-0, 50  DO  BEGIN 

READF, 1,R2 

OI  (J)-R2 (3) 
ENDFOR 
CLOSE, 1 
RETURN 
END 
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PRO  GETVE,  ALT,  WK,  VLBH, VSECPOS , DENOM 

GETXRATE ,  ALT,  A1P ,  A3S ,  B3P  ,  C3P ,  W3D 

GETDENS, OALT, ODEN 

A=.5 

B00=.125 

QO=3E-ll 

Q05=1.69E-1 

0=INTERPOL (ODEN, OALT, ALT) 

VLBH=A1P 

VSECPOS=C3P/2. 

DENOM=(l.+QO*0/A) 

WK=  (  (A3S+B3P+C3P/2  .  )  *B00*Q05)  /DENOM 

RETURN 

END 

0GETXRATE 

GGETDENS 


THIS  PROCEDURE  USES  GETXRATE. PRO  AND 
GETODEN.PRO  TO  GET  VOLUME  EXCITATION  RATES 
FOR  VK,  LBH,  AND  SECPOS  FOR  SLNTINT.PRO 
WRITTEN  BY  K.L.  MCELROY 
15  MAY  1992 
MODIFIED: 14  SEP  92 


;THIS  IS  1/LIFETIME  OF  A  STATE (S"-l) 

;THIS  THE  BRANCHING  RATIO 
;THIS  IS  O  QUENCHING  COEFF . (CMA3  SA-1) 
;THIS  IS  THE  FRANK-CONDON  FACTOR  FOR  0-5 
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;THIS    PROCEDURE    CONVERTS    VOL    EX 
;FROM    GETVE.PRO    TO    SLANT    INTEN 
;FOR   VK,LBH,    AND    SECPOS    OF    N2 
/WRITTEN    BY    K.L.MCELROY 
;23    APR    1992 
/MODIFIED:    15    MAY    92 


pro    slntint,attO,alt,svk,slbh, ssecpos, denom 

get ve, alt, vvk, vlbh, vsecpos, denom 

att-alt*0+attO 

delalt-alt (0)-alt(5) 

tophvk-ab3 (delalt/alog (vvk  (5) /vvk (0) ) ) 

topvevk-vvk (0) *exp (-2 . /tophvk) 

tophlbh-abs (delalt/alog (vlbh  (5) /vlbh (0) ) ) 

topvelbh-vlbh (0) *exp  (-2 . /tophlbh) 

tophc3pi-abs (delalt/alog (vsecpos (5) /vsecpos (0) ) ) 

topvec3pi-vsecpos (0) *exp (-2 . /tophc3pi) 

volslnt, vvk, alt, at t, topvevk, tophvk, svk 

volslnt, vlbh, alt, at t, topvelbh, tophlbh, slbh 

volslnt, vsecpos, alt, att , topvec3pi, tophc3pi,  ssecpos 

svk=svk*le-6 

slbh-slbh*le-6 

ssecpos«ssecpos*le-6 

return 

end 

Ggetve 
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;THIS  PROCEDURE  CALCULATES  THE  SLANT 
/INTENSITIES  FOR  01(2972). 
; WRITTEN  BY  K.L.  MCELROY 
;11  JULY  1992 
/MODIFIED:  16  SEP  92 

PPO  OISLNT, ATTO, ALT183, SOI 

ATT-FLTARR(18  3) 

ATT-ATT*0+ATTO 

A2972-0.045  /THESE  ARE  THE  EINSTEIN  COEFFIENTS 

A5577-1.06 

B-.l  /THIS  IS  THE  BRANCHING  RATIO 

K2-2.1E-8 

K3-5E-11 

V-'  ' 

READ, 'WHAT  VERSION  OF  FILES  DO  YOU  WANT?',V 

/THIS  SECTION  GETS  THE  FIRST 
/PRODUCTION  FACTOR  PI 

OPENR, 1, ' 90089FLX'+V+' .PRT' 

STR1-'  ' 

WHILE  STRPOS (STR1, 'VOLUME  PRODUCTION  RATE')  EQ  -1  DO  READF,1,STR1 

WHILE  STRPOS (STR1, 'SPECIES   01')  EQ  -1  DO  READF,1,STR1 

FOR  J-1,3  DO  READF,1,STR1 

ALT-FLTARR(51) 

01S-FLTARR(51) 

Rl-FLTARR(ll) 

FOR  J-0, 50  DO  BEGIN 
READF, 1,R1 
ALT(J)-R1  (0) 
01S(J)-R1  (3) 

ENDFOR 

CLOSE,  1 

/THIS  SECTION  GETS  THE 
/DENSITIES  FOR  02+  &  ELEC 

OPENR, 1, ' 90089ION' +V+' .DAT' 

A-ASSOC(l,FLTARR(183)  ) 

02P-FLTARR(183) 

ELEC-FLTARR(183) 

32P=REVERSE (A (5) )  /MAKES  02P  AND  ELEC  ARRAYS  GO 

ELEC=REVERSE(A(6) )  /FROM  HIGH  TO  LOW 

CLOSE,  1 

/THIS  SECTION  GETS  N4S  DENSITY 

OPENR, 1, ' 90089NOX' +V+' .DAT' 

A=ASSOC(l,FLTARR(183) ) 

N4S-FLTARR(183) 

N4S-REVERSE (A (1) )  /MAKES  ARRAY  GO  FROM  HIGH  TO  LOW 

CLOSE,  1 

/THIS  CONVERTS  01S  FROM 

/A  FLTARR(51)  TO  FLTARR(183) 


77 


ALT18  3-FINDGEN(183)  *2  .+8  6 

ALT183-REVERSE(ALT183) 

01SDEN«INTERPOL(01S,ALT,ALT183) 


GETDENS, OALT, ODEN, 02DEN, N2DEN, TEMP 
ODEN-REVERSE (ODEN) 
02DEN-REVERSE (02DEN) 
N2DEN-REVERSE (N2DEN) 
M-ODEN+02DEN+N2DEN 


;THIS  SECTION  GETS  0,02  6  N2 
/DENSITIES  AND  TEMP  NEEDED  FOR 
;P4 

; THESE  COMMANDS  MAKE  THE 
/DENSITIES  GO  FROM  HIGH  TO  LOW 


T-REVERSE(TEMP) 
K4-4.7E-33*  (300/T) A2 


P1-01SDEN 

P2-K2*B*02P*ELEC 

P3«K3*02P*N4S 

P4-K4*ODENA2*M 

L-A2972+A5577 

VOI«(Pl+P2+P3+P4)/L*A2  972 


DELALT-ALT183 (0) -ALT183 (5) 

TOPHOI-ABS (DELALT/ALOG (VOI (5) /VOI (0) ) ) 

TOPVEOI-VOI (0) *EXP (-2./TOPHOI) 

VOLSLNT, VOI,ALT183,ATT,TOPVEOI,TOPHOI, SOI 

SOI-SOI*lE-6 

RETURN 

END 

8GETDENS 


; MAKES  TEMP  GO  FROM  HIGH  TO  LO 

;THIS  SECTION  CALCULATES  THE  V 
;FOR  012972 


;THIS  SECTION  CALCULATES  THE 
; SLANT  INTENSITY  FOR  012972 


;THIS  CONVERTS  TO  RAYLEIGHS 
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